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> ■ Abstract 
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■ We discuss the interpretation of Euclidean correlation functions at finite tem- 

. perature (T) and their relationship with the corresponding real-time Green's 

0^ 



functions. The soluble 2+1 dimensional Gross-Neveu model in the large-A^ 
^ limit is used throughout as a working example. First, the real-time bound 

Q_i' state, identified as an elementary excitation at finite T, is solved. The bound 

state mass, the dispersion relation at low momenta, the coupling constant 
D . and decay constant are calculated. To characterize the structure of the bound 

' state the on-shell form factor is carefully introduced and calculated. Then we 

examine the corresponding screening state and contrast the screening mass, 
^ ' coupling constant, decay constant and the screening Bethe-Salpeter ampli- 

' tude with the real-time quantities. We find that, although they can be used 

as qualitative indicators in the low-T regime, the screening states at finite T 
in general do not reflect the properties of the corresponding real-time bound 
states. Besides, other relevant issues, such as the subtlety of the real-time 
manifestation of conservation laws due to some internal symmetries at T 7^ 0, 
the temperature dependence of the pseudoscalar spectral function and its sum 
rule, and the high-T limit of the screening state and its implications to the 
dimensional reduction, are also discussed in detail. 
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I. INTRODUCTION 



The formulation of QCD, or any other field theory, at finite temperature is well- 
established and several equivalent approaches exist. Nevertheless, since finite temperature 
systems are determined not only by their ground states but also by all the excited states, 
the study of finite temperature field theory really involves different physics with respect to 
the zero temperature case and different concepts need to be introduced 

In the existing literature the real-time dynamics is mainly discussed in a perturbative 
context. However, we are often interested in non-perturbative real-time physics. Typical 
non-perturbative methods for QCD calculation at finite T are the operator product expan- 
sion and lattice simulations. Both these approaches are intrinsically Euclidean and limited to 
static properties. The possibility of extracting real-time information from these Euclidean 
quantities needs to be addressed. Due to the lack of experimental data, we do not have 
enough intuition to guide our approaches. Therefore, it is particularly useful to have some 
exact results which could provide us with insights into such matter. However, demanding 
exact solutions often implies a sacrifice in the direct phenomenological relevance. 

In this paper, we address these issues in the U(l) chirally symmetric Gross-Neveu 
model in 2 + 1 dimensions and in the large- iV limit: this soluble model allows us to 
perform exact calculations both for real-time and screening dynamics. This model is often 
used in the literature whenever exact results are required and, besides, some relevance 
to QCD is also desired. In fact, the Gross-Neveu model has several features in common 
with QCD. For instance, this model undergoes a chiral phase transition at finite T, as it is 
expected to happen in QCD, providing a good place for a qualitative modeling of the temper- 
ature dependence of real-time dynamics. However, some caution is needed: the Gross-Neveu 
model is not confining and, therefore, not every physical interpretation of the singularities 
associated to quark and antiquark can be literarily applied to QCD. 

More specifically, we solve the real-time bound state and study its properties as functions 
of temperature. We calculate its mass, its coupling constant, its decay constant and its 
spectral function. In addition, we carefully characterize the size of this bound state by 
calculating its on-shell form factor. Similarly, we also solve the corresponding screening 
state by calculating the screening mass, coupling, decay constant, and the screening Bethe- 
Salpeter amplitude. The choice of screening observables, and of the corresponding real-time 
ones, has been suggested from the available lattice calculations . Eventually, we contrast 
the real-time observables and the corresponding screening ones. 

While some of the results presented here can be found in the existing literature, and we 
include them only to make the presentation of our points coherent and self-contained, others 
are less well-known. Among these new, or insufficiently discussed, topics we would like to 
emphasize the following. 

(1) Due to the lack of Lorentz invariance at finite T a general amplitude A{u!,p) has often 
very different functional dependences on u"^ and p^. An obvious consequence of this non- 
uniformness is, for example, the non-covariant energy-momentum dispersion relation for 
elementary excitations. A less obvious but more interesting consequence is related to how 
conservation laws due to some internal symmetries manifest themselves in various physical 
processes. The Goldberger-Treiman relation and the effective charge of the pion are used as 
illustrations. 



2 



(2) To characterize the structure of a bound state at finite T in a physical way we intro- 
duce the on-shell form factor in the elastic limit. The peculiarity of the relevant kinematic 
condition and the existence of an additional form factor are carefully addressed. Then we 
isolate the singularity structure that characterizes the on-shell form factor as a function of 
the spatial momentum transfer and, therefore, determine the size of the bound state. 

(3) We derive an exact sum rule, which states that the zeroth moment of the pseudoscalar 
spectral function is temperature independent. 

(4) The screening Bethe-Salpeter amplitudes in scalar and pseudoscalar channels are calcu- 
lated explicitly. 

(5) We obtain the asymptotic formula for the screening mass in the high-T limit and show 
that this result demonstrates the failure of dimensional reduction in this model. The reason 
of this failure is also explained. 

Strictly speaking, the 2+1 dimensional Gross-Neveu model with a continuous symmetry 
has no finite-T phase transition beyond the leading order in 1/iV, due to the severe infrared 
singularity associated with massless Goldstone bosons. One might question the relevance of 
our results at leading order in This problem can be avoided by working in 2 + e spatial 
dimensions (0 < e < 1), where the infrared catastrophe is absent. Since the limit e — > is 
smooth at the leading order in our results presented in this paper should be understood 
as such. We could have also addressed the same questions in the 3+1 dimensional Nambu- 
Jona-Lasinio model [0. However, since our present interest is more related to answering 
conceptual questions rather than finding direct phenomenological applications, we prefer to 
use the 2+1 Gross-Neveu model because of its renormalizability. Moreover, we believe that 
the qualitative physics associated with the chiral restoration at finite T is similar in all these 
models, independently of their dimensionality. In any case, we shall take notice of those 
results that are directly linked to the specific dimensionality. 

The paper is organized as follows. First, we give, in Sec. 0, a brief pedagogical review 
of how real-time dynamics is formulated at finite T. Then, in Sec. [m| , we introduce the 
model and collect all those results that are useful in the rest of the paper. The real-time 
pseudoscalar bound state and its properties are calculated in Sec. In this same section we 
also discuss the Goldberger-Treiman relation, the on-shell form factor and derive a sum rule 
for the spectral function. The screening state is discussed in Sec. 0, where we also contrast 
minutely the real-time bound state and the screening state. The last subsection of Sec. |V| is 
dedicated to the high temperature limit of the screening state and to the implications that 
this result has for dimensional reduction. We reserve Sec. [V| to a summary of our work and 
to conclusions. 



II. GENERAL INTERPRETATION 

In this section we review the connection between correlation functions, which are the 
typical output of theoretical calculations, and the real-time response of a system to external 
perturbations. The following material is most likely well-known to experts [Q; nonetheless, 
since we are specifically concerned with the proper interpretation of Euclidean correlation 
functions, we find useful to include it here as a convenient reminder. 

An experimentalist perturbs a system with an external probe, and then measures the 
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effect of this perturbation with some detector. In a given theory we describe this external 
perturbation using an interaction Hamiltonian of the form 



^ext(t) = / dxJ{i,{t, a;))Kxt(t, x) , with Kxtlt < 0, a;) = , (1) 

where Vey±{t,x) is the external perturbation, which is switched on at t = 0, and couples 
to the system through the current J{ilj(t,x)) (the hat indicates operators). The system is 
described by a state, or ensemble of states of the Hamiltonian, and the measurement of the 
response to the perturbation is performed through another current, which can be the same as 
the one that couples to the perturbation. For example, we may perturb the system with an 
external electric field, which couples to the charge density, and measure the resulting change 
in the same charge density. Therefore, we typically measure the change of the expectation 
value of a given current due to the presence of the perturbation: 

s{j{m))^{JiHme.t-{-mmo, (2) 

where {A)ext is the expectation value of the operator A in the state (or ensemble of states) 
describing the system in the presence of the external perturbation, and (^4)0 the correspond- 
ing expectation value in the unperturbed system. Here and in the following we drop the 
spatial label for simplicity. 

If the perturbation is weak, we can expand in the perturbation, and keep only the term 
linear in the external field. Then the linear response of the system to the perturbation is 
simply proportional to the retarded correlation function: 

S{J{m)) = / dt'G^it - t')Kxt(t'), (3a) 
J —00 

where the retarded correlation is defined as 

^ ,(,) T.(e-^../n.(.,W) .W,(0))1) ^^^^ 

Tr (e'^o/T\ 

and iplt) = e^^^°'i[j{0)e~^'^^° , Hq is the unperturbed Hamiltonian, and we consider a system 
in thermal equilibrium at temperature T. The physical content of this equation is known as 
the fluctuation-dissipation theorem, i.e. correlation functions of a system describe not only 
the correlations (fluctuations) of the system in a given channel, but also its linear response 
(dissipation) to weak external perturbations. 

We say that the system has a (real-time) elementary excitation in the channel described 
by a given current, if we get a resonant (very large) response to the corresponding external 
perturbation at a given frequency. If we rewrite Eq. (|3aD in frequency-momentum space 
(restoring the spatial dependence) 

5(J)(^,fc) = G«(a;,fc)V;xt(^,fc), (4) 

we can see that the appearance of a pole in the correlator as a function of uj gives this kind 
of resonant response, which will be called the bound state. It is clear from Eq. (11) that the 
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resonance frequency depends on spatial momentum k: cu — u!{k). The function a;(fe) is the 
so-called dispersion relation for the excitation. 

Alternatively, we can use a time-independent external perturbation (a; = 0), and study 
the static response of the system. In this case, we are interested in the modification of the 
system in spatial distribution, e.g. how the effect of the perturbation dies out with distance 
(screening). In momentum space there can still exist purely imaginary poles in the retarded 
correlator, but now they arc interpreted as pure exponential terms in the decaying response. 
At zero temperature only, due to Euclidean invariance for rclativistic theories, there is a 
one to one correspondence between these poles in the imaginary momentum and the ones 
in frequency. 

The resonant frequencies of the system (poles of the propagator) are important charac- 
teristics of the system, but they are not the only information we get from experiment (and 
from the propagator). The strength with which the external probe couples to the system is 
also important. For instance, excitations of the system that do not couple to the external 
probe may not be relevant. In the propagator, the information about the coupling strength 
is carried by the residue of the pole at the resonant frequency. 

A nice example of the concepts we have just reviewed is the classic electron plasma 
probed through the charge density: J{t,x) = ip^{t,x)ip{t,x). The response of the plasma 
to a static point charge is just the well-known Debye screening 

S{J{x)) (X Jdt G^{t, x) oc e-"^^l="l , (5) 

and the corresponding pole in the static propagator. 



1 

is the Debye mass: m|, = STrrioe^/ kT {no and e are the electron density and charge, respec- 
tively). 

But the plasma also possess real-time excitations, plasmons, with characteristic frequency 
ujl^{k = 0) = 47re^no/me, where nie is the electron mass. These excitations are described, 
near the plasmon pole, by the Green's function 

G^iuj, k) « ^— — - , (7) 

u - u;ph(fc) + i-f 

where 7 is the damping rate which gives a finite width to the plasmon. The frequency- 
momentum (or dispersion) relation for plasmons is uJp^i^) = '^ph(O) + 

From the preceding example it is quite clear that in general real-time excitations and 
screening states describe different physics. Nevertheless, there is a connection between 
measurements in Minkowski and Euclidean space. In general, causality and unitarity make 
possible the following dispersion relation for the correlation function 



G 



R 



{u,,k)='-rdu , (8) 

TT J -00 u — ui — le 



where the spectral function p is proportional to the total cross-section in the given channel. 
It is important to realize that the dispersion relation at finite T is in frequency only and the 
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spatial momentum is treated as a parameter. In the plasmon example, if the plasmon were 
the only excitation of the system, the spectral density p would just be the on-shell condition 
for the plasmon: 

p{u,k) = \^e{u)5{u'-ul^{k)). (9) 

The power of the dispersion relation is that the physical content of the correlation is entirely 
embodied in the spectral function p, which determines the correlation in the whole complex 
plane, apart from a possible additive polynomial in uj. For instance, if we analytically 
continue the correlation function to imaginary frequencies, we obtain 

~g^{u,k) = - r du^^^. (10) 

TT J-oo U — iuj 

This analytically continued correlation function calculated at the Matsubara frequencies 
ujn = {2n — l)7rT (for fermions) is the discrete Fourier transform of the correlator calculated 
in the Euclidean region (finite temperature correlator): 

Tr fe-^<'/^J(V'(r))J(V^(0))]) 
g{r) = / . ' ''^^ (11) 



where ijj{T) = e'^^°i/j{0)e~^^'-' . This formula is valid in the interval < r < 1/T, while 
for T outside this interval we can use the periodicity condition, which for fermions reads 
g{T + 1/T) = —g{T) (it follows from the periodicity of the trace and the anticommutation 
property of the fermionic fields). 

Vice versa, if we know the correlation function at the discrete set of Euclidean points 
specified by the Matsubara frequencies, the analytic continuation u;„ — >• —iu ± e provides 
us with the retarded/advanced correlation for continuous real-time frequency. This con- 
tinuation is unique only when the expression of the correlation function is explicitly free 
of singularity in the entire cj-plane except on the real axis and obeys certain convergence 
property at Itul ^ oo p|. 

Finally, due to the lack of Lorentz invariance at finite T, a general Green's function 
G{uj, p) often has very different functional dependence on uj"^ and on p^. For example, G{0, 0) 
strongly depends on the order in which we take the limits uj —>■ and p 0, i.e. G{0^, 0) = 
lim^^Qlimp^Q G{uj,p) is usually different from G{0,0~^) = limp^olim^^Q G(^u,p). (The 
innocent looking function [x + y)/{x — y) in the limit oi x, y —>■ illustrates the relevant 
subtlety.) Of course, the physics of the specific process selects which of the two order is 
relevant. In particular, G{0'^, 0) is intrinsically Minkowskian, because it represents a process 
whose characteristic time is much shorter than the heat-bath response-time. We shall then 
call this kind of processes fast processes. Correspondingly, G{0, 0"*") is intrinsically Euclidean 
and the related processes are called slow processes, since the thermal environment has enough 
time to respond to the external perturbation. As we will see later, conservation laws due to 
internal symmetries only hold in fast processes. In slow processes the apparent violation of 
conservation laws should be understood in the sense of the grand canonical ensemble, not 
as the violation of the fundamental physics laws. 
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III. THE GROSS-NEVEU MODEL 



In this section we introduce the Gross-Neveu model in 2 + 1 dimensions, and present 
its main features at finite temperature. Corresponding formulae for the model in 1 + 1 
dimensions can be found in Ref. 



A. Lagrangian 



The Gross-Neveu model is defined by the Lagrangian density 



L = tpi'y ■ dip + 



2N 



(12a) 



where ip is & 4-component Dirac spinor with the color indices implicit. Equivalently, we can 
write 



C = ipij ■ dip — ip(a + m'^^)ip — — :t(o"^ + tt 

^9 



(12b) 



where a and vr are the auxiliary scalar and pseudoscalar bosonic fields, respectively. We 
study this model in the limit N —>■ (yo with the coupling constant fixed. This Lagrangian 
is invariant under a continuous U{1) chiral symmetry that is dynamically broken at low T. 

B. Effective potential and mass generation 



The standard approach of the effective potential at finite temperature yields the 
critical temperature and the dynamically generated mass. To the leading order in 1/N, the 
effective potential for a and tt fields is given by the one-loop expression 



TV 



2N 



(2^ 



+ 00 



(13) 



where ujn = (2n — l)7rT. Up to an irrelevant constant the sum over n can be done yielding 

Kff(cr,0) = ^a2-2iV j ^{v/pT^+2rin[l + exp(-v/pT^/T)]} . (14) 

In this last and the following equations we have used the symmetry of the Lagrangian, 
and rotated the (a, vr) field in the a direction making vr = 0. Alternatively, we can read a 
in the following formulae as representing a/o^ + tt^. The integral is divergent, and a high 
momentum cutoff A is required. We make V^g finite by adding a counterterm of the form 



-CT = -ir(o- + TT ) 

2 V TT TT 



(15) 



i.e. we have defined the bare coupling constant 'n / g\ = ^/S'^ ~ + ^5 with g"^ = 5'^(k) the 
renormalized coupling constant. 
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Then, the renormahzed effective potential is 



TT I 3' ' 2 



dz In 



1 + exp(-v^) 



(16) 



where fi = k — 7r/g'^{K,), which is independent of the renormahzation scale k, will turn out 
to be the dynamical quark mass at T = 0. 

The quantum theory is determined by the fluctuations around the minimum of the ef- 
fective potential, which we find from the stationary condition 



dVes{cr,0) 



da 



N 

—CTm - + 2T lu [1 + exp(-crm/T)] } 







(17) 



When T > Tc (high T phase), the expression in braces is positive definite, and the only 
solution to Eq. ([TtD is am = (symmetric phase). The critical is the value of T for which 
the expression in braces is zero with am = 0, i.e. = /i/ In 4. 

When T < T^, the expression in braces is zero for the following value of am'- 



am{T) = /i 1 









+ 2T In 







1 -4exp(-/i/T) 



(18) 



which is the absolute minimum of the effective potential, Eq. ([T6|), while the other stationary 
point, (T = 0, is a maximum for these values of T. Consequently, the symmetry is dynamically 
broken and the corresponding quark mass is m(T) = am(T). 

In Fig. |1| we plot the dynamically generated mass as a function of temperature, i.e. 
Eq. (p^Sl): it goes exponentially to fi in the zero temperature limit, while it approaches zero 



as 2^ ii{Tc — T) for T approaching Tc from below. 

For later convenience we define the following function of T, which is proportional to the 
first derivative of V(,s{a, vr) with respect to a"^ evaluated at (a = m, tt = 0), 



R{T) = fi- m{T) - 2T ln[l + exp(-m(T)/T)] 



1 



if T < Tc 

ifr>Tc 



(19) 



C. Pseudoscalar Bubble Graph 

In this model, the only one-particle irreducible loop graph (in the large N limit) that is 
relevant for our purpose is the pseudoscalar bubble graph Hp. In order to extract the bound 
state and screening states, we need this graph in Minkowski space for general [uj,p). We 
first use the imaginary-time formalism to perform the loop integral: 

iYip{u,p)=NT / T^Tr ^75T ^1^T , (20) 

n=^^J (27r)^ ( k--f-m k ■ - p ■ - m j 

where k = {iujn, k), p = {iu!,p), ujn = (2n — l)7r with n = 0, ±1, ±2, ■ ■ -, and u = 2ttI with 
/ = 0, ±1, ±2, ■ ■ -. The sum over frequencies can be performed by using standard contour 
integral techniques |TD[, and we obtain: 
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/-A d^k tanh(Efc/2T) / 4k ■ p 
Iip{u,p) = -Nj ^ — 



Ak ■ p 



[p 



p p^ + 2k ■ p 



(21) 



ko=Ek 



where E^. = + . The above equation is infinite and can be made finite by adding the 
contribution from the counterterm zIIct = d'^Lc-^jdo^ = —N{A — n/g'^ — fi)/^. We then 
analytically continue u into the entire complex plane. 

For calculating the bound state mass, it is sufficient to have the expression for zero 
external momentum in the region < tu^ < 4m^ for the pseudoscalar channel: 



- + -R{T) + 



9^ 



TC 



An 



E 

dE tanh — - 



1 



2T E^-uj^/A 



- + -{ R{T) + - In 



"2m + ui' 




2m — uj. 


A Jm 



1 + E^-uj^/A^ 



(22) 



In the region u"^ > Am?, which is needed for the spectral function, the bubble develops 
an imaginary part, and the corresponding formulae are: 



Re [ilip{u;)] 



+ - <! R{T) + - In 
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vr 



UJ + 2m 
UJ — 2m. 



UJ^ 

-—R dE 



A Jm 1 + e^/^ E^ - uyA ^ 



lm[iUp{uj}\ = — a;tanh — 



(23a) 
(23b) 



where Ff indicates the principal part of the integral. 

For studying the screening phenomena we need the pseudoscalar bubble, Eq. (pOD, in the 
Euclidean region for uj = 0, i.e. in the static limit. It is convenient to expand the hyperbolic 
tangent by using the formula 



tanhx 



E 



x2 + [n{n - 1/2)Y 
and perform the fc-integral. The resulting expression is 



^lUpituj ^0,p) = ^ + R(T) -T Y: 

9 n=— o> 

where MI=ujI + m^{T). 



\ AMI + 



In 



AMI +P^ + Vp^ 



AM^+p^-^ 



(24) 



(25) 



D. Two-point function 



To the leading order in the retarded two-point function for pseudoscalar current 
involves only a sum over a geometric series of the pseudoscalar bubble graph: 



W^>)t{^.P) = I dxdt exp{tujt-tp-x)e{t){[Mt,x),M0,O)])T 



iUp{uj,p) 
itgyN)Up{uj,p) ' 
(26) 
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where J^{t, x) = 'ilj{t, x) 275 ilj{t, x), and we have renormahzed the two-point function in such 
a way that its lowest order coincides with the renormahzed bubble. 

When the above equation has a pole in real frequency a; at a given spatial momentum 
p, this pole is identified as the elementary excitation or the bound state. In the static limit 
(u; = 0), on the other hand, the lowest singularity for —p^, denoted by m^, corresponds to 
the screening mass. 

IV. REAL-TIME DYNAMICS 

In this section the real-time dynamics is studied. We first solve the bound state in the 
pseudoscalar channel, including the bound state mass, with its dispersion relation for small 
momenta, the coupling and decay constants. The Goldberger-Treiman relation is examined 
carefully. We then calculate the on-shell form factor and the charge radius. Finally, the 
spectral function and its sum rule are derived. 

A. Pseudoscalar Bound State 

As discussed in section |I|, a physical excitation in a given channel is signaled by a pole 
in the corresponding correlation function. Therefore, a bound state in the pseudoscalar 
channel (there is no elementary pseudoscalar excitation in the Lagrangian) exists, if there 
is a solution to the equation obtained by equating to zero the denominator of Eq. ( p6D 

= 1- j^illp{uj = m^) 

\2m + m^] ml p 2 1 

= i?(T) + -ln[^--— J--y^ ^^TT^E^3^' (27) 

where we have used the expression for the pseudoscalar bubble, Eq. (p^, valid for < < 
4m^. This equation has always the solution m,r = 0, which is lower than the energy of the 
unbound particles 2m, as long as T < Tc (-R(T) = 0). The study of the residue associated 
with this solution can also help in determining the fate of the bound state when T — > T^. 

As long as we are only interested in the bound state mass, the bubble calculated at zero 
momentum p is sufficient. However, it is also instructive to calculate the energy as a function 
of momentum. In fact the system at finite temperature has a preferred reference frame, the 
heat bath, and we expect an explicit loss of covariance for the energy-momentum dispersion 
relation. To this end, we perform the calculation of the bound state energy retaining terms 
up to the first order in p^ and a;^ (a;^ is also small in the limit of small p^). The bubble in 
this limit can be written as: 

where we have dropped terms of order p^, u^, p^uj"^ and higher. It should be emphasized 
that when obtaining Eq. (^) the bubble is expanded in powers of p^ first, and then the 
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corresponding coefficients are expanded in powers of uJ^ . Again by solving the bound state 
equation ill = N/ g"^, we find 



uj\p')=v\T)p' + 0{p'') (29a) 



where the function 



2 ^ ^^dx{x ^ + 3x ^)tanh(ma:/2r) 
^ ^ ' 2j^dxx-Hanh{mx/2T) ' ^ ^ 

has the following limits: lim^-^o '^^^(^) = 1 , as it should, and limy^T;, f^(T) = 1/2. The fact 
that f^(T) < 1 implies that the speed of the pion is reduced by the thermal environment, 
even though the pion is still massless. 

As stated in section ||, the bound state coupling constant to a quark-antiquark pair is 
given by the residue of the full pion propagator on the mass-shell 

which is also proportional to the residue of the pole of the correlation function, Eq. (^), at 
u = m^. 

The dependence of g1{T) on temperature is shown in Fig. |^. When T 0, m —>■ fi, and 
gl — > 4:Hfi/N (zero temperature limit). When T — >• T~, m(T) — > 0, and gliT) 0"*". In this 
last case, the coupling gl{T) approaches zero logarithmically: lim^^^,- gl{T) oc (logm)~^ oc 
[log (Tc — T)]~^. In contrast, the corresponding coupling in the four dimensional Nambu- 
Jona-Lasinio model in the chiral limit approaches a finite constant at — 0~^, and then it 
jumps to zero for T > Tc 0. 

The bound state solution in the pseudoscalar channel for T = has the same energy of 
two free quarks, and one might wonder about the fate of the bound state in the limit T — > 
T~. The coupling constant result gives us a clear answer suggesting that the bound state 
disappears at the phase transition. More precisely, we should conclude that the pseudoscalar 
meson decouples from its constituents at the phase transition. 



B. The Goldberger-Treiman relation 

We can also introduce a pion decay constant defined as the residue of the axial current 
correlation function at the pion pole (this definition at T = yields the usual expres- 
sion (0|J^|7r(p)) = —ip^f-w). So we need to calculate the axial-pseudoscalar bubble graph 
A^{iuj,p), defined by replacing in Eq.(^) one of the 275 with 7/^75. It is easy to verify that 
A^{iuj,p) =Pf,A{iuj,p) with 

^ ,. NT^ f\ m , ^, 

Anuj,p) = > / da— -. — r :— , (31) 

where uj = 2ttI is the bosonic Matsubara frequency. The above equation can not be naively 
used to analytically continue iuj into the entire complex plane, since iu has singularities off 
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the real axis. However, for calculating the pion decay constant we only need A at p = 0. 
Then the a-integral can be done, while the sum over n can be carried out by the standard 
contour integral technique, yielding, 

, N , / X \ m 
A{iuj, p = 0) = — dx tanh — — , (32) 



2nJrr. \2TJ x2-(itu)V4 

which now can be trivially continued and then put on mass-shell. Hence, the pion decay 
constant is given by 



MT) = g^{T) A(0+, 0) = g^{T) — —tanh j , (33) 

which is valid in the broken phase. Since there is no pion state in the symmetric phase, f-n{T) 
is not defined when T > T^. The vanishing of fj^ at T = merely refiects the fact that the 
chiral symmetry is restored. Using the explicit expression of g'^{T), the above equation can 
be rewritten as 

g^{T) UiT) = 2m(T) , (34) 

which is the Goldberger-Treiman relation at finite T. Notice that gA = ^ identically in the 
large- limit. 

If we were not careful in using A(0+, 0) to define fT^iT) we would have got an expression 
of fn{T) not satisfying the Goldberger-Treiman relation. For example, if we put the axial- 
pseudoscalar bubble graph on the mass-shell before taking the limit p — *^ 0, we would have 
effectively used A{0,0^) = ^tanh(m/2T) in Eq. (|33|) . The physical reason that the pion 
decay should be regarded as a fast process is that the decay happens instantaneously and 
the thermal environment does not have enough time to respond. 

At this point it is appropriate to comment on the fact that it is not trivial how, at 
finite temperature, the conservation laws that are the consequence of symmetries of a the- 
ory manifest themselves in the real-time dynamics. Ward identities associated with global 
symmetries can be straightforwardly generalized from the T = case to finite T in the 
Euclidean formalism. However, it is very subtle how to subsequently analytically continue 
these relations between Euclidean Green's functions to corresponding relations between real- 
time quantities and at which stage one should put the external lines on their mass-shells. 
Therefore, one can not simply assume that Euclidean Ward identities immediately apply 



also to the on-shell quantities, as was done in Ref. |TT|]. It appears that the validity of the 
conservation laws depends on the specific mathematical prescriptions, but this dependence 
is not arbitrary: it refiects the nature of the formulation of field theories at finite T in terms 
of ensemble averages. Only in fast processes, where the heat bath does not have time to 
respond, conservation laws hold explicitly. Whereas in slow processes, where a measurement 
always involves the feedback of the heat bath, conservation laws are no longer manifest. We 
will encounter an example of such a "violation" in the next subsection. 

C. On-shell Form Factor 

At any temperature the on-shell form factor in the elastic limit provides a physical 
measure of the size of an elementary excitation. In fact, even at finite T, we can imagine 
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the following experiment. In the heat-bath frame, we scatter a lepton off the pion. Just as 
in the usual zero temperature scattering, we select those leptons that have scattered off an 
on-shell pion rather than off something else by using the appropriate kinematic conditions. 
In our case, this is possible because the pion is massless, while an unbound quark- ant iquark 
is massive. The elasticity here is important, since the kinematics selects out slow processes 
in which the scattered lepton sees not only the pion but also the thermal cloud around it. 
The size defined this way provides a snap shot of the elementary excitation immersed in an 
equilibrated thermal environment. 

The on-shell form factor has also been studied in Ref. []T2|. However, the author appar- 
ently did not realize several complications that arise when defining on-shell form factors at 
finite T, such as the proper on-shell condition at finite T and the appearance of additional 
form factors. 

On general ground, a typical three-point function with two pion lines and one photon 
line, as shown in Fig. ^ has the following structure in the heat-bath frame 



The additional form factor ■ p'] vanishes identically at T = and it is related to 

the heat bath. To minimize the environmental effect and hence to emphasize the intrinsic 
structure of the pion it is natural to use the the invariant function F to define its size; to 
this end, we can select the spatial components fi{p,p'), which receive contribution only from 
F, i.e. we consider a "magnetic scattering process". 

For low energy elastic scattering, the momentum transfer is q = p' — p = {0,p' — p), 
i.e. po = Pq and = p'"^. Due to the breaking of the explicit Lorentz invariance at finite 
temperature, the on-shell condition for the pion does not imply p"^ = uj'^ip) — p^ = fn^. 
Therefore, there exist two independent variables in the on-shell form factors F and G, even 
in the elastic limit. We choose these two independent variables to be uj{p) and = {p' —pY- 
Then the radius of the pion at finite temperature is defined as 



The factor of 4, rather than the usual 6, is due to the fact that we are in two spatial 
dimensions. This definition coincides with the usual definition at zero temperature. The 
overall normalization of F[qr^, C(j(p)], which depends on T because the charge of the pion 
gets screened {F[q^ = 0,u{p)] < 1 in general), does not affect the charge distribution. The 
limit of uj{p) m^r is to ensure that the radius is measured in the heat-bath frame. This 
condition need not be imposed at zero temperature, due to the Lorentz covariance. In case 
one has trouble to imagine a charge radius for a massless particle, we could have included 
a small quark mass in the Lagrangian such that > 0. This formal process is not really 
necessary since the limit — is smooth. 

The graph that gives the leading contribution to F[q'^, uj{p)] is shown in Fig. The 
advantage of probing the pion by a photon instead of a scalar, the a, is that meson-pole 
dominance is avoided, consistent with our intention to emphasize the intrinsic structure. 
Using the pion-quark-antiquark coupling derived in the last section, we have 



/^(P, P') = {Pf. + p'^,)F[p\ p'\ p-p'] + 6^^G[p\ p'\ p-p']. 



(35) 




(36) 
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Mp,p')=9i{T)T y: 



(Pk Tr75(/c ■ 7 + m)'y^{k ■ 7 — J9 ■ 7 + m)'-fi{k ■ 7 — p' ■ 7 + m) 



(2vr) 



[k'^ — m?] [{k — pY ~ "^^1 [(^ " P'Y ~ "^^1 



(37) 

{iun, k), 



where m = m{T) is the mass gap. In the elastic hmit the momenta are k 
P = (Po,p) andp' = {po,p') with po = uj{p) = v{T)\p\ + 0{\p\^). 

Since we are deahng with a slow process, we should enforce the mass-shell condition at 
the very beginning. We can take advantage of the facts that the pion is massless and that 
the scattering is elastic and enforce the on-shell condition by setting all external Matsubara 
frequencies to zero. This step is perfectly legal, because the Euclidean point iu = i27il 
with I = happens to lie on the real axis and coincides with the point a; = 0. Then, the 
calculation of the on-shell form factor from Eq. (^7|) involves the following steps: (1) combine 
the denominators using the Feynman parameter representation; (2) carry out the spatial 
momentum integral; (3) eliminate all momentum variables in terms of the two independent 
variables and uj{p); (4) extract F[q'^,u{p)] from fi{p,p') using Eq. (P^D; (5) take the limit 
uj{p)— > m^. Without going through the details, which are tedious but straightforward, we 
directly give the final result 



F[q',uip)^0] 



9l{T) 



27r 



T 



00 



da 



1 



(38) 



The Matsubara frequency sum can be readily carried out, yielding 



F[q\u{p)-.Q]=gl{T) 



tanh(m/2T) 
A-Krn 



{m/T) 



sinh (m/T) 



0{q' 



(39) 



Using the explicit expression of g1{T) in Eq. (|30|) we immediately find that the effective 
charge F[Q,uj{p)] < 1, with the equality sign only at T = 0. The fact that F[q^ ,uj{p)] 
vanishes in the limit T ^ Tc merely reflects the fact that the cross section for a lepton to 
hit the pion also vanishes in the same limit, consistent with the decoupling of the pion from 
its constituents at the chiral restoration point. As another consistency check, F[0,C(j(p)] can 
also be calculated from the full pion propagator S'^[c<j,p^] through the Ward identity, in the 
hmit of — s> 



F[QMp)] = 9l{T) 



d 



dp 



,S-'[uJ,p'] 



(40) 



U)=Ll)(p) 



which is a consequence of the residual static gauge invariance at finite temperature. In the 
above equation one must use the full pion propagator in the limit of slow processes, i.e. 
taking a; = before sending p^ to zero as it was done in Eq. (pSf); this limit differs from 
the one used to produce Eq. (|28D which is defined for fast processes. As we anticipated, the 
charge conservation that guarantees F[0] = 1 at T = no longer holds at finite T, due to the 
well-known charge screening. The effective charge F[0] is the pion charge plus the charge 
in its thermal cloud. The explicit charge conservation only can be expected in a deeply 
inelastic process, in which the on-shell condition is enforced at the end of the form factor 
calculation. Of course, the total charges of the test particle and the heat bath is definitely 
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conserved. However, in a slow scattering process, the lepton only sees the charges localized 
near the test particle, but not those charges far away at the spatial boundary. 

It would be interesting to calculate directly the thermal-cloud charge distribution induced 
by a pion. However, this calculation appears to be beyond the scope of the linear response 
theory. The reason is that, according to the linear response theory, the induced thermal 
cloud is proportional (in momentum space) to the charge distribution of the test particle, 
and the temperature dependence appears only in the proportionality function, i.e. the 
retarded correlation function. In our case, the test particle itself has a non-trivial internal 
structure, which is subject to change under the thermal environment. Therefore, there is 
back-reaction between the induced thermal cloud and the structure of the pion, which is not 
an external probe. This intricate entangling of the test particle and the thermal environment 
necessarily prompted us to the calculation of the on-shell form factor. 

Now the radius of the pion in the heat-bath frame can be obtained easily 



The complete curve of {r'^)T as a function of T is displayed in Fig. ^, which makes sense only 
in the symmetry broken phase where the pion exists. At zero temperature (r^)o = l/(3/i^), 
which can also be verified directly from the covariant calculation at T = 0, where F[q^\ = 
[4m^/(— g^)]^/^ sin~-^[(— g^)/(4m^ — q"^)]^^"^, with = ql — q^. Near critical temperature 
(r^)r l/(18^c) = 2(ln2)^/(9/i^). Although the pion's size decreases as temperature 
increases from zero to the critical point by about factor of 3, the characteristic size scale 
of the pion remains to be This is related to the fact that the pion, once it forms, is 

always a tightly bound state when T < T^. 

It is also interesting to point out that, in general, the radius {r'^)T, defined through slow 
processes, can not be simply estimated by the threshold of the triangle graph, which is a 
standard practice at T = 0. In fact, if one did so, one would obtain the wrong estimate 
(r^)r ~ which diverges when T —>■ Tc. The phase space integral, which produces the 

factor in the brackets in Eq. (^), can play an important role at finite temperature. As can 
be seen clearly from Eq. (^), the actual singularity structure of F[q'^,uj{p) — > 0] in is 
always controlled by the lowest threshold 4M^ = 4(7r^T^ + m^), which stays of order /x^ in 
the entire broken phase. A physical interpretation can be given for Eqs. (^) and (^T]): the 
thermal environment tends to wash out coherence beyond the thermal wave length (1/T) in 
slow processes. 

At this point the physical picture is clear. At low temperature, the pion state is very 
much the usual pion state with exponentially small correction from the thermal environ- 
ment. As temperature increases the quark and antiquark become less and less likely to bind 
together to form the elementary excitation. However, once they bind together, the elemen- 
tary excitation remains qualitatively similar to its zero temperature counterpart. Above 
the critical temperature the pion completely resolves into its constituents and no longer 
exists. So the phase structure of the system is strongly reflected from its spectral content in 
Minkowski space. 

In sharp contrast, the screening Bethe-Salpeter amplitude defined in the Euclidean space 
does not share this feature, since it is insensitive to the real-time singularity structure, as we 




1 r 



1 - 



m/T 



(41) 



sinh(m/T)- 
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shall demonstrate in the next section. Nonetheless, the screening Bethe-Salpeter amplitude 
indeed provides a qualitative measure of the size of the corresponding bound state, though 
only in the phase where the bound state can be identified as an elementary excitation. This 
is due to the fact that the on-shell form factor and the screening Bethe-Salpeter amplitude 
have thresholds proportional to . 

Can the result in Eq (^) be an accident of the 2+1 dimensional Gross- Neveu model? 
To settle this question we need to study the on-shell form factor in 3+1 dimensions, i.e. 
in the Nambu-Jona-Lasinio model. In 3+1 dimensions, we find that the resulting formula 
is not very different from the one in 2+1 dimensions, Eq. apart from some trivial 

factors, which are independent of the momentum transfer, the only modification is that the 
integrand of the a-integral is raised to the power 1/2, as it is obvious from dimensional 
counting. Therefore, the lowest singularity in q^-plane is still given by 4Mf Although 
now the form factor requires a cutoff in the Matsubara frequency sum, the charge radius 
is independent of this cutoff. This implies that the charge radius in 3+1 dimensions has 
qualitatively the same behavior as in the 2+1 dimensional case. Our result, a pion charge 
radius that remains finite at the critical temperature both in the 2+1 Gross-Neveu model 
and in the 3+1 Nambu-Jona-Lasinio model, is in sharp contradiction with the result found 
in Ref. [|T^], a pion charge radius that diverges at the critical temperature in the Nambu- 
Jona-Lasinio model. Since the result of Ref. [0 is not expressed through a simple formula 
such as our Eqs. ( ^8]) and (|4T|) , but is obtained by means of a numerical integration of an 
expression that involves principal value definitions, we have not been able to pin down the 
source of this discrepancy. However, it is worth mentioning that Schulze used the temporal 
component of the triangle graph to define the charge radius, which we believe is less desirable 
due to the contamination from the second form factor G associated with the heat bath. It 
is also not clear at what exact stage he enforced the on-shell condition, and hence it is not 
clear whether his form factor is defined for fast or slow processes. Either of these two points 
could lead to very different results and, more importantly, to different interpretations in 
terms of physical measurements. 



D. Spectral Function 

We now calculate the spectral function in the pseudoscalar channel. We give the result 
for p = 0, where the cut contribution near the origin, due to the scattering of the current and 
thermal particles, vanishes. Below Tc we have a contribution from the continuum {uj > 2m) , 
and a pole contribution from the bound state: 



TT 



9' 



6{uj' - mi) + e{uj' - 4m'{T)) peont(o^) 



(42) 



where the bound state mass m^, = 0, and the coupling gl(T) is given in Eq. (0). Above 
Tc only the contribution from the continuum survives. The continuum part of the spectral 
function is related to the real and imaginary parts of the bubble graph, see Eqs. (0), in the 
following way 

/ \ Im[inp(c<j)] 

^ (1 - i9yN)Re mpH]f + ii9VN)lm [iUpH]?' 
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We show the spectral function at different temperatures in Fig. ^. The most remarkable 
feature is that near the phase transition the shape of the spectral function is very different 
from its shape at low temperature. A large peak develops just above threshold. The peak 
actually diverges right at the phase transition, though its integrated strength stays finite. At 
the same time, as expected, the pole becomes weaker and the threshold goes to zero. These 
qualitative features are not artifacts of the particular model we are considering here. They 
actually are the characteristics of a continuous phase transition and critical phenomena. For 
example, the strong peak near the threshold is a reflection of the fact that the corresponding 
susceptibility diverges at T^. These characteristics would still persist to a certain extent when 
the phase transition is weekly first order. Furthermore, as we argued in Ref. |T^, QCD shares 
all these qualitative features. 

The strong peak near the origin in the spectral function, when T is close to Tc, indicates 
a resonance- like excitation. This resonance corresponds to a complex pole of the two-point 
function on the second Riemann sheet of the cj^-plane, with a very small imaginary part. 
Using explicit formulae in Sec. |T|, one can analytically continue the two-point function onto 
the second Riemann sheet. Then, it is not hard to find that the real and imaginary parts of 
the pole in the symmetric phase are proportional to t/ ln(t~^) and —t/ ln^(t~^), respectively, 
where t = [T — T^/Tc <C 1. As T gradually increases from this resonance pole moves 
almost parallel to the real axis from the origin to the right initially; and it eventually marches 
into the first quadrant. When T is outside the scaling region, the imaginary part of this 
complex pole becomes so big such that it does not makes sense anymore to call it a resonance 
excitation. 

It appears that the diminishing of the pion strength and the magnification of the quark- 
antiquark continuum strength in the spectral function near the critical region provides a 
possible physical picture for the failure of the standard a-model scenario in predicting 
the nature of the chiral-symmetry-restoration phase transition in the Gross-Neveu model 
in 2+1 dimensions. To substantiate this statement a detailed investigation is necessary. The 
reason is that the a field in the relevant cx-model does not necessarily correspond to any 
real-time bound state pole, but more likely to the resonance pole mentioned above. 



E. Sum Rule 

To quantitatively characterize the weakening of the pole and the growing of the peak 
above threshold we derive an exact sum rule for the temperature dependent part of the 
zeroth moment of the spectral function. 

To derive the sum rule we need to find out the asymptotic behaviors of the two-point 
function in the deep Euclidean region and of the spectral function in the large-o;^ limit. 



Using the expression in Eq. (|2ll) it is easy to obtain the pseudoscalar bubble graph in the 
Euclidean region 

.np(..Q.p ^ 0) ^ ^ H- ^{fi(r) - fun- §- 4/;,^^^^} . (44) 

and its asymptotic form at large-Q^ 
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^.n.w,.^o,^-f.(l.^)-±^^.o,«-^ (45) 

where ((■ ■ ■)) stands for the thermal average 

/ dE-—^A{E). (46) 

Jm L -\- C I 

For example, {{E"^)) = 3({3)T^, in the symmetric phase where m = 0. The corresponding 
two-point function, defined in Eq. (EB), is 



JM ( 

(J5J5)tW,0) 



N r 4 



Afi 16/^2 64/i3-167r2(m3 + 3((E2)))/3 



+ OiQ-") . (47) 

On the other hand, the spectral function has the following asymptotic form at large u;^, 
as can be derived from Eqs. ( p3D and 



,^^)..^^U.^^-mVUo(.-^). (48) 

To study the temperature dependence we make the subtraction, following Ref. |T3], 

A(J5J5)(^Q,0) = {J,J^)t{iQ,Q) - {J,J,>)t'{iQ,^) = r du^ y^t^a ' (^9) 

with Ap(c<j) = pTiyj) — Pt'{^)- Since the leading term in pxiyj) is independent of T, as seen 
from Eq. (^8]), we immediately have Ap ~ u;~^. This implies that the zeroth moment of Ap 
exists. Using Eq. (|3), and hence 

A(^.^.>(.Q.0)^g^ ^""^y«^^>» -.O(Q--). (50) 
we arrive at the following expression appropriate to derive the sum rule 

where A(m3 + 3((E2))) = (m^ + ?,{{E'^)))^ - (m^ + 3((E2)))^,. Taking the limit ^ 00 we 
get the desired sum rule 

roo 

/ rfa;^ Ap = , (52) 

Interchanging the order of the dispersion integral and the ^ 00 limit is allowed since the 
zeroth moment is finite. 

We have proved that the very same sum rule also holds in the 1+1 dimensional Gross- 
Neveu model and QCD [|13|, though we have used a totally different derivation in those 
cases and the corresponding sum rules would only converge logarithmically. It appears that 
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the temperature dependence of the zeroth moment of the spectral function often shows 
quahtatively similar behaviors, no matter what is the underlying microscopic physics. 

Another interesting point, which was especially clear in the derivation we used for the 
1+1 dimensional Gross- Neveu model and QCD [T^, is that the zeroth moment of the spectral 
function could be related, in asymptotically free theories, to expectation values of appro- 
priate operators via the operator product expansion. That result cannot be immediately 
applied to the 2+1 Gross- Neveu model, since this model is not asymptotically free and 
possesses a nontrivial ultra-violet fixed point. Therefore, it would be extremely interesting 
to investigate whether a formalism similar to the operator product expansion |]TB[ can be 
extended to this case. 



V. SCREENING PHENOMENA 

As discussed in section screening phenomena are associated with responses to time- 
independent external perturbations. Therefore, we only need to consider static Green's 
functions and purely spatial correlation functions. Static Green's functions are essentially 
Euclidean and hence, in contrast to the time-dependent responses we have discussed in the 
preceding section, well-suited for the lattice formulation. In this section we study the static 
two-point and three-point functions that relate to screening masses and screening Bethe- 
Salpeter amplitudes, respectively, since these quantities are studied on the lattice [^|[. 



A. Screening mass 



The pseudoscalar screening mass m is the lowest solution to the equation 



—Ilp{iuj 0,ip) 



^\r{T) + 2T V 

TT ' 



tan 



-1 



(53) 



~ 2 

mi. 



In the high temperature 



When R{T) = (T < TJ, the solution is obvious: p'^ = 
phase R(T) ^ 0, the solution can only be obtained numerically. 

Contrary to the naive expectation that higher than lowest Matsubara modes decouple 
at high T, we point out that the frequency sum in Eq. ( |53D cannot be truncated without 
introducing an error in the screening mass that is not of order 1/T. This fact suggests that 
dimensional reduction does not take place in this model [^]. Details on this point will be 
presented later. 

We can define a coupling constant for the screening state analogous to the one defined 
for the real-time bound state case 



' d_ 
dp- 



Ilp{iuj 0,p) 



(54) 



Then, near the screening-mass pole, the static pseudoscalar propagator behaves like 
g1{T)/{p^ + m^). In Fig. Q gliT) is plotted in dashed line. Even though g1 agrees with 
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gl at T = 0, due to the Euclidean invariance, it becomes drastically different from gl as 
temperature increases. In particular, does not vanish in the symmetric phase, but rather 
grows linearly with T in the high-T limit. In the broken phase gl{T) has a particularly 
simple form: gl{T) = Anm/[N tanh(m/2T)]. 

As in the real-time case, an analogous "decay constant" for the "on-shell" screening pion 
can be defined. Again it is easy to verify 



UT) = UT) ^ tanh 



(55) 



which also vanishes in the symmetric phase. It is interesting to notice that the "Goldberger- 
Treiman" relation holds exactly in the screening case 



UT)UT) = 2m{T) 



(56) 



We believe that this result is not accidental, because the screening state in the static limit 
can be regarded as a bound state in some 1+1 dimensional theory at zero temperature with 
the same symmetry as the original theory. This Euclidean "Goldberger-Treiman" relation 
is in agreement with the results of Ref. |]ll] . 

In order to explicitly see the manifestation of the chiral symmetry also in the screening 
masses, it is instructive to calculate the scalar screening mass. The scalar bubble graph (in 
the static limit) is easily related to the pseudoscalar counterpart 



NT 



TT 



E 



4m^ 



p2 \ 4M2+p2 



In 



4M2 + p2 + ^ 



4M2 + p2 _ ^ 



(57) 



Similarly to the pseudoscalar screening mass, the scalar screening mass is found by solving 
the equation 







3_ 

N 



Us{iuj 0,ip) 



-^\r{T) + 2T V 



p2 — 477^2 



p2 ^ AMI -P^ \ 4M2 - p2 J • 



(58) 



When T < Tc, R{T) = and the solution is again easily obtained p^ = Am? = 777.2 (notice 
that = is not a solution in this channel). For T > Tc the dynamical quark mass 
w?(T) vanishes, and the scalar screening mass equation becomes identical to the one for the 
pseudoscalar screening mass. 

In Fig. ^ we display numerical results for the screening masses in both channels, pseu- 
doscalar and scalar. In the low temperature phase, chiral symmetry is broken, and the pion 
is the relevant Goldstone boson: this fact is also reflected in the screening mass. In the high 
temperature phase, chiral symmetry is restored, and we find the expected parity doubling, 
i.e. the degeneracy of the screening masses in the parity mirrored channels. 

We already know that the screening state is very different from the real-time bound 
state. One might, however, naturally ask whether the screening state be somehow related 
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to the resonance excitation near Tc- In particular, one might want to introduce some non- 
trivial external momentum dependence in the screening state and hope that, by making 
this external spatial momentum complex, the screening pole might be continued into the 
resonance pole. While this option is in principle possible, we think that a more relevant 
question is whether the knowledge of the static screening state as a function of temperature 
is sufficient by itself, i.e. when we do not have complete knowledge of the analytic structure 
of the two-point function, to provide enough direct information on the resonance state. We 
think that the answer to this latter question is unlikely to be affirmative. The reason is that 
the functional dependences of the two-point function on frequency and spatial momemtum 
at finite T are genuinely different, due to the explicit lack of Lorentz invariance. Therefore, 
it is not possible to probe, in general, the full information carried by the frequency variable 
by varying a spatial component of the momentum. In addition, the following three points 
also support this general view. First, the screening pole is purely imaginary for any T, 
whereas the resonance pole is complex in general. Second, the screening state is well defined 
for any T, whereas the resonance is only identifiable near the critical region. The third point 
is more specific to the model we have consider: when T ^ Tc, the screening mass squared is 
proportional to t, with t = {T — Tc)/Tc ^ 1, whereas the resonance pole is proportional to 
t with logarithmic corrections (see section IV. D). 



B. Screening Bethe-Salpeter amplitude 



Another observable that is calculated on the lattice at finite temperature is the screen- 
ing Bethe-Salpeter amplitude (SBSA). Our definition of this amplitude in the pseudoscalar 
channel is 



(j){x, y) = dr dz {ij{T, x, z)-i^i){T, x,z + y) ?/'(0, 0, 0)75^(0, 0,0)). 



(59) 

This definition differs from the usual lattice definition by an extra integration over r. We 
perform this extra integration for computational convenience, since it allows us to consider 
only the zero frequency mode. In terms of physical content of the SBSA, this difference 
should be irrelevant. To the leading order in the Feynman graphs contributing to the 
SBSA are those shown in Fig. ^ Then the Fourier transform of 0(r, y) has the following 
form 



0(Pl,P2) 



<^l(Pl,P2) 



l-gmp{iuj^Q,p={p^,Q))/N 
where 0i(pi,p2) is the off-shell screening BS amplitude 
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0i(pi,P2) = iNTY^ 



d^k 
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An explicit calculation yields 
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Therefore, the complete SBSA is 



" y 27r ' 7 271^ l-^72np(p = (z^-.0,pi,0))/iV ^^^^ 

Typically, we are interested in the dependence of the SBSA on y for large x, so that we 
project on the lowest "screening state". In the limit x ^ 00 the pi-integral picks out the 
lowest singularity, p\ = < 4:M^, and the SBSA becomes 



The scale that controls the size of the screening wave function is clearly y Mf — / 4, as 
one can show by deforming the integration contour to the imaginary axis. Since Mf = 
(ttT)^ + m^, and = at low temperature, in the low T limit this scale is m oc /x. At high 



temperature m = 0, and rriTr oc T, therefore in this limit the scale is set by y tt^T^ — m^/4 ~ 
T. 

As we have argued in the preceding section, a better measure of the spatial distribution 
of the quark and antiquark inside a meson is provided by its on-shell form factor. The SBSA 
qualitatively measures the size of the elementary excitation, or bound state, in the entire 
broken phase. However, the SBSA above the phase transition does not yield any information 
about the nature of the corresponding (nonexistent) real-time bound state. 

We also calculate the SBSA in the scalar channel using exactly the same procedure we 
have used in the pseudoscalar channel, and we find very similar results. The relevant mass 



scale is now given by yMf — rri^/4, which has the same high temperature limit as the 
scale for the pseudoscalar channel. In the broken phase (T < Tc), since = 2m, and 
Ml = (vrT)^ + m^, ^ — rfi^ / 4 = ttT. This last result is consistent with the fact that 
there is no binding energy in the scalar channel, and the only screening is due to the thermal 
mass. 



C. Real-time vs. screening phenomena 

Now we have all the ingredients to make an explicit comparison between the real-time 
bound state and the screening state. 

Mass: When T < (the chiral symmetry is spontaneously broken), both the real-time 
and the screening pion masses vanish, rriTriT) = rh.„(T) = 0, while the scalar mass fhaiT) = 
2m{T). However, when T > Tc (chirally symmetric phase), the real-time pion decouples and 
the spectral function has contributions only from the quark-antiquark continuum, whereas 
the screening pion mass m.,t{T) becomes degenerate with the screening scalar mass ificr(T) 
and both of them asymptotically grow linearly with T (see Fig. 

Coupling: The real-time coupling of the pion to the quark-antiquark state, gl(T), and 
the corresponding screening coupling, gl{T), have very different behaviors as functions of 
temperature, as it is clearly shown in Fig. |^, and they only coincide at T = 0. In particular, 
the screening coupling gl(T) does not provide any signal for the fact that the real-time 
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excitation decouples when T > Tc, in fact gl(T) does not show any distinctive feature near 
the critical region. 

Decay constant Both /7r(T) and fn(T) start at the same value at T = and gradually 
decrease to zero at T = Tc. The vanishing of fniT) and fniT) when T > Tc is due to the fact 
that the axial current decouples from the pseudoscalar current when the chiral sjnumetry is 
restored. This constraint forces a similar behavior for both f-n{T) and f-KiT), even if their 
numerical values between T = and T = Tc are different. Above the phase transition, /^r 
loses its meaning, and fn{T) remains zero. 

Size: As seen from Fig. ^, the sizes defined through the on-shell form factor in the elastic 
limit, (r^)^ and the screening Bethe-Salpeter amplitude, 1/Mf, have qualitatively similar 
behaviors in the broken phase. The reason for this similarity is that they share the same 
singularity structure in low spatial momentum transfer. Nevertheless, we find again that, 
while loses its meaning above Tc, the size defined through the screening Bethe-Salpeter 
amplitude does not give us any signal of this disappearance of the real-time state. 

We believe that most of the above-mentioned qualitative features will also be present in 
3+1 dimensions, although some of the details presented here are certainly specific to this 
model. Therefore, we can safely conclude that the screening observables, in general, do not 
necessarily reflect the behavior of the corresponding real-time observables. These two kinds 
of observables are often qualitatively different. 



D. High-T limit and dimensional reduction 

The present model study can also offer us very interesting information about the possi- 
bility of dimensional reduction for screening Green functions at high temperature. 
In the symmetric phase, equations ( |53D and (|^) coincide and can be rewritten as 

OO 

;^ - ^ = 4 E tan-i , (65) 

Tc T ^1 J4{2n - 1)2 - x2 JA{2n - 1)2 - 



where x'^ = p'^/{t[TY. In the high-T limit we can solve the above equation iteratively by 
setting X = Y^'^^Ci (/i/vrT)*. In particular, the first two coefficients can be easily found by 
solving 

^ = 4 E .... 'V. . tan- ^ , (66) 



Tc t^i JA{2n - 1)2 - d JA{2n - 1)2 - d 



and 



4-[ Co 4(2n-l)2 Co ^ 

ci vr;i'44(2n-l)2-cr [4(2n-l)2-c2]3/2 J4(2n - 1)2 - eg ' ' ^ ^ 



Therefore, the asymptotic form for the screening masses in the high-T limit is, 

^7r,o- = Co ttT + ci /i + 0(Ai/T) . (68) 
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Numerically, cq ~ 0.982 and ci ~ —0.947. Eq. ( |68D is plotted in dashed line in Fig. y. 
A salient feature of this asymptotic behavior is that the high temperature limit of the 
screening mass is Cq ttT, where the numerical value of Cq has contributions from all the 
Matsubara frequency. If dimensional reduction took place in the high-T limit, we would 
find that only the lowest (or at most a finite number) modes would contribute to cq. In fact, 
any truncation in the frequency sum in Eq. (^), and then in Eq. (|66D, would give a different 
Cq: the resulting error in the screening mass would not be suppressed in the high-T limit. 

The reason for the failure of the dimensional reduction picture is that the 2+1 dimensional 
Gross-Neveu model lacks the necessary scale hierarchy: there is no elementary bosonic 
particle whose zero mode could dominate over other modes of order T, and the coupling 
constant cannot provide a second lower scale compared to T, as it happens in asymptotically 
free theory []T6[, since this model has a nontrivial ultraviolet fixed point. In fact, the high- 
momentum behavior of the dimensionless coupling G{k,) = ng'^in) has the form /i/k), 
with G* = TT (the value of G* is scheme dependent). This nontrivial fixed point implies that 
this dimensionless coupling "runs", in the high temperature limit, to a finite value (since 
K ~ T |T^). Therefore the effective interaction strength at high T between all the Matsubara 
modes (~ TG{T), once the fields have been appropriately re-scaled) is strong and its strength 
in units of T becomes independent of temperature (~ TG*). 

In this regard the Gross-Neveu model in 1+1 dimensions is very different from the one 
in 2+1 dimensions. The 1+1 dimensional model is asymptotically free and g'^{T) runs to 
zero logarithmically, providing the scale hierarchy between T and g'^(T)T ||16[. In this model 
dimensional reduction takes place and screening quantities are well reproduced by the lowest 
modes that are weakly interacting: the leading contribution to the screening mass is given 
by the free theory, cq = 2, the first correction to this non-interacting behavior is of the order 
g^(T) and it is correctly reproduced by the reduced theory. Since g^(T) runs to zero the 
reduced theory becomes exact asymptotically. 

Finally, we should point out that the concept of dimensional reduction that we have 
discussed in this subsection and in Refs. [p!6| , p!7| , is different from the concept of dimensional 
reduction employed for predicting the nature of phase transitions at finite T in the a-model 
scenario |T^. Our criterion for dimensional reduction is more stringent and applies to a 
different temperature regime. We require that the static Green's functions of the original 
theory at high-T limit be matched by the corresponding Green's functions in the reduced 
theory to a given accuracy, whereas in the cr-model scenario one requires that the reduced 
theory be able to match the thermodynamical singularities near the critical region and hence 
give the same critical exponents. Therefore, the success/failure of one of the two dimensional 
reduction pictures does not directly implies the success/failure of the other. 



VI. CONCLUSION 

We have shown in a specific example, the Gross-Neveu model in 2+1 dimensions, that 
the direct connection between the real-time bound state and the corresponding screening 
state is lost at finite temperature in general, especially near and above phase transition. 

In particular, the real-time pion disappears after the phase transition, as it is signaled by 
the vanishing of its bounding energy and of the coupling to its constituents g'^{T). In spite of 
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this, the corresponding screening state is still strongly "bound" in the high temperature limit, 
and the screening coupling g1{T) grows linearly with T. The screening mass approaches the 
asymptotic value of 0.9827rT, which is much smaller than the sum of the masses of two 
noninteracting fermions, 2'kT. 

Moreover, the screening Bethe-Salpeter amplitude yields sizes of order 1/T both in the 
scalar and pseudoscalar channels at high T, when the real-time bound state no longer exists. 
Therefore, no relevant information about the nature of the real-time pion can be inferred 
from this screening amplitude. On the contrary, we can give a natural definition of the size 
of the real-time pion (below the phase transition) by means of the on-shell form factor in the 
elastic limit, since the on-shell form factor has direct information of the relevant coupling 
gl{T). Furthermore, we show that the associated size cannot be estimated solely from simple 
threshold considerations, reflecting the fact that the thermal environment tends to wash out 
coherence beyond thermal wavelength in slow (compared to the equilibration time of the 
thermal bath) processes. 

We have also computed the exact pseudoscalar spectral function of this model. The most 
important features of this spectral function are the diminishing of the strength of the pion 
pole and the magnification of the strength of the quark- ant iquark continuum near the phase 
transition. In addition, we have derived an exact sum rule for this spectral function: the 
integrated strength of the spectral function is independent of T. 

Many of the features of the 2+1 dimensional Gross-Neveu model that we have listed 
above are dominated by the chiral symmetry and should have their analogues in the real 
world. For instance, very similar results were found in the 1+1 dimensional Gross-Neveu 
model |T^, and in the Nambu-Jona-Lasinio model [0 in 3+1 dimensions. On the other 



hand, certain specific aspects certainly depend on the details of the theory. For example, we 
find that the 2+1 dimensional Gross-Neveu model does not undergo dimensional reduction 
in the high-T limit. But the possibility of using a lower dimensional theory to describe 
static correlations of the original theory when T is large is strictly tied to the short distance 
behavior, which is very different in this model from the one in QCD |jl7 . 



Finally, we also find that it is instructive to classify physical processes according to their 
characteristic time scales relative to the time scale required by the thermal environment 
to respond to an external perturbation. We have explicitly illustrated, in the cases of the 
Goldberger-Treiman relation and of the effective charge of the pion, that conservation laws 
due to internal symmetries are manifest only in fast processes but not in slow ones. 

It is our hope that all the lessons learned from this exactly soluble model study will 
provide some useful insights to QCD. 
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FIGURES 



FIG. 1. The dynamical mass as a function of temperature. The mass and T are in units of /x, 
which is the dynamical mass at T = 0. 

FIG. 2. The bound state coupling constant, g^(T), as a function of temperature (solid line). 
The dashed line is the coupling constant for the screening state g'^{T). Both couplings are in units 
of their zero temperature value: Ann/N. 

FIG. 3. The spectral function in the pseudoscalar channel at fe = as a function oiuj/ ji: (a) in 
the broken phase at T = (dashed line) and T = 0.9524 Tc (solid line); and (b) in the symmetric 
phase at r = 1.25 Tc (solid line) and T = 2.5 Tc (dashed line). The arrows in (a) denote the pion 
poles at cj = 0, with a slight displacement for visual clarity. The height of the arrows indicates the 
relative strengths of the poles. 

FIG. 4. Feynman graph for the three-point function, whose on-shcU value gives the form factor. 
The thick external lines are the pion lines and the wiggle line is the photon line. 

FIG. 5. On-shell charge radius, normalized to its value at T = 0, {r'^)T/{r'^)o, as a function 
of temperature in the symmetry broken phase. For comparison, the "size" estimated from the 
screening Bethe-Salpeter amplitude, 1/Mf , is also plotted in dashed line. 

FIG. 6. Screening masses in the scalar (cr) and pseudoscalar (tt) channels in units of /x. It is 
also plotted in dashed line the asymptotic formula at high-T limit, rhcr = rhT^ = 0.982 ttT — 0.947/x. 

FIG. 7. Feynman graphs for the screening Bethe-Salpeter amplitude. 
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